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Abstract 

We analyze the Higgs intense-coupling regime, in which all Higgs particles of the 
Minimal Super symmetric Standard Model are light with masses of the same order and 
the value of tan f3 the ratio of vacuum expectation values of the two Higgs fields is large, 
in the framework of Supergravity scenarios with non-universal soft Supersymmetry 
breaking scalar masses in the Higgs sector. In particular, we calculate the relic density 
abundance of the lightest neutralino candidate for cold dark matter and the rates in 
direct and indirect detection at present and future experiments. We first show that 
while in the mSUGRA model this regime is disfavored by present data, there are regions 
in the parameter space of models with non-universal Higgs masses where it can occur. 
We then show that because of the large value of tan f3 and the relatively low values 
of the neutral Higgs boson masses, the cross section for neutralino-nucleon scattering 
is strongly enhanced in this regime and would allow for the observation of a signal 
in direct detection experiments such as CDMS-Soudan. The expected sensitivity of 
gamma-ray detectors like GLAST might be also sufficient to observe the annihilation 
of neutralinos in such a regime. 
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1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) [1], the scalar sector is extended to 
include two Higgs doublets fields to achieve the breaking of the SU(2)l x U(1)y electroweak 
gauge symmetry. This leads to the existence of five Higgs particles: two CP-even Higgs 
bosons h and if, a CP-odd or pseudoscalar Higgs boson A, and two charged Higgs particles 
[1, 2]. At the tree-level, the spectrum is determined by two basic parameters: the mass 
of the pseudoscalar Higgs boson Ma and the ratio of the vacuum expectation values of the 
Higgs fields, tan/9 = V2/V1. If the pseudoscalar Higgs particle is very heavy. Ma ^ Mz-, one 
is in the so-called decoupling regime [3] in which the lightest CP-even h particle is SM-like 
and has a mass that is close to the Z boson mass while the other CP-even Higgs particle H 
and the charged bosons are very heavy and degenerate in mass with the pseudoscalar A 
boson, Mh ~ Mh± ~ Ma- The MSSM Higgs sector reduces then to the one of the SM, but 
with a Higgs particle that it rather light: M^^^ ^ Mz at the tree-level but M™^"" ~ 100- 
140 GeV, depending on tan (3 and the strength of the important radiative corrections which 
need to be included [2, 4]. In the opposite scenario, in which Ma < Mz [5] called the anti- 
decoupling regime in Ref. [2], it is the lighter CP-even h which is degenerate in mass with 
the A boson, while the heavier H boson is SM-like with a mass close to M^™ ~ M™^. 

An intermediate scenario, the so-called intense-coupling regime [6], is characterized by 
a rather large value of tan (3 and a mass for the pseudoscalar A boson that is close but 
not equal to the maximal (minimal) value of the CP-even h (H) boson mass. In such a 
scenario, an almost mass degeneracy of the neutral Higgs particles of the model occurs, 
Mh ^ Ma Mh ~ 100-140 GeV, while the couphngs of both the CP-even h and H 
particles to gauge bosons and isospin up-type fermions are suppressed, and their couplings 
to down-type fermions, and in particular 6-quarks and r leptons, are strongly enhanced. 
The interactions of both Higgs particles therefore approach those of the pseudoscalar Higgs 
boson which does not couple to massive gauge bosons as a result of CP invariance, and for 
which the couplings to isospin — ^ (+|) fermions are (inversely) proportional to tan/3. 

This scenario leads to a very interesting collider phenomenology which has been discussed 
in detail in Ref. [6]. In particular, it has been shown that one has to face a rather difficult 
situation for the detection of these particles at the LHC, since the branching ratios of the 
usual interesting decays which allow the detection of the CP-even Higgs bosons are too small 
and cannot be used anymore, while the dominant bb and t'^t~ decay modes have too large 
backgrounds. Even when using some rare Higgs decays, the detection of the three individual 
Higgs bosons is very challenging in general and in some cases even impossible at the LHC. 
At the future International Linear e+e" Collider (ILC), the three peaks can be resolved but 
the measurement of the masses of the particles is more difficult than in other scenarios. 

The studies on the intense coupling regime mentioned above have been performed only in 
the framework of an effective low-energy MSSM in which the parameters which break softly 
Supersymmetry (SUSY) are incorporated by hand, leading to a Higgs sector that is practi- 
cally disconnected from the SUSY particle sector. However, it would be more interesting to 
consider constrained SUSY models, which are theoretically more appealing and which have 
a smaller number of basic input parameters. A scenario that is widely used as a benchmark 
for constrained Grand Unified SUSY Theories (SUSY GUTs), is the minimal Supergravity 
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model (mSUGRA) [7] in which there are only five input parameters: a universal value mo for 
the soft SUSY-breaking masses of all the scalars, a universal soft SUSY-breaking gaugino 
mass parameter M1/2, a common trilinear Higgs-sfermion coupling [the three of which are 
defined at the GUT scale, Mqut ~ 2 x 10^^ GeV], tan (3 and the sign of the Higgs-higgsino 
mass parameter /i. An interesting question is whether the intense coupling regime can occur 
in such a scenario. 

A second question that one might ask is whether this scenario is compatible with cos- 
mological observations and, in particular, with the observed amount of dark matter in the 
universe. Indeed, it is now well established that luminous matter makes up only a small 
fraction of the observed mass in the universe and a weakly interacting massive particle, 
identified in the MSSM with the stable lightest ncutralino, is one of the leading candidates 
for the "dark" component which is needed to explain astrophysical data on the rotation 
curves of large scale structures of the universe [8] and which has been recently measured via 
the anisotropies of the cosmic microwave background by the WMAP satellite with a very 
good precision [9]. A related question is whether in such a scenario, the lightest neutralino 
of the MSSM would lead to observable effects in the present and near future experiments 
which are devoted to the direct and indirect detection of these dark matter particles. 

In this paper, we attempt to answer to these questions. We show that in the framework 
of the mSUGRA model, the particle spectrum is so constrained that the intense coupling 
regime is realized only in a very limited area of the parameter space; in addition, because 
it occurs only for extremely large values of tan /3, the obtained spectrum could potentially 
lead to a conflict with experimental data from collider physics. We then investigate the 
Non Universal Higgs Mass (NUHM) model, discussed in Refs. [10, 11] for instance, in which 
the equahty of the soft SUSY-breaking masses of the two Higgs fields m^i and mj^j with 
the common mass term mo for squarks and sleptons at the GUT scale is relaxed. Such a 
non-universality structure might, for instance, occur in SUSY GUTs in which the fields Hi 
and H2 belong to different multiplets [10]; it can also be obtained in the low-energy limit 
of some phenomenologically appealing string scenarios such as heterotic orbifold models 
[12] and effective heterotic Anomaly Mediated SUSY Breaking scenarios [13]. The NUHM 
model is then equivalent to an effective low energy MSSM where, in addition to the mSUGRA 
parameters, on has the higgsino mass parameter and the pseudoscalar Higgs boson mass 
Ma as free input parameters. This additional freedom will allow to realize the intense Higgs 
coupling regime, while complying with all known accelerator and astrophysical data, in a 
much larger region of the parameter space compared to the mSUGRA case. 

The paper is organized as follows. In section 2 we review the NUHM model and discuss 
the impact of non-universal scalar mass parameters on the Higgs sector as well as on the 
cosmological relic density of the dark matter lightest neutralino and on its detection in direct 
and indirect searches. In section 3, we analyze the Higgs intense regime in the NUHM model: 
we first explain the procedure to obtain the spectrum and the various constraints that we 
impose on it, and then present two concrete examples in which this regime is realized; we 
then discuss the implications of these scenarios for the relic density of the neutralino and 
its direct and indirect detection rates. For completeness, we briefly discuss in section 4, the 
case of the mSUGRA model. A brief conclusion is given in section 5. 
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2. The NUHM model 



2.1 The impact of non— universality on the Higgs sector 

Let us first recall that in the mSUGRA model, in addition to the four basic continuous input 
parameters, the common soft SUSY-brcaking scalar mo and gaugino M1/2 mass parameters 
and trilinear coupling Aq, defined at the GUT scale, and the ratio of vcvs tan/3, the sign of 
the Higgs-higgsino parameter /j, is also free. This is due to the fact that /x^ is determined by 
the minimization of the scalar Higgs potential which leads to electroweak symmetry breaking 
(EWSB). At the tree level, one has in terms of tan /3 and the two soft SUSY-breaking Higgs 
mass terms defined at the SUSY breaking scale Ms — 0{1 TeV), 

^ " tan2/3-l 2^ ' 

which, for reasonably large values of tan/3 ( > 5), can be approximated by 



~ -ml^ - -Ml , (2) 



1 

2^ 

with being negative to effectively generate EWSB. Again at the tree level, the mass of 
the CP-odd Higgs boson is approximately given by 



= V , (3) 

which can be rewritten using eq. (2) 

Ml ^ mjj^ - mjj^ - Ml , (4) 

The masses of the neutral CP-even Higgs bosons and the charged Higgs particle are then 
obtained from the usual tree-level relations 



Ml + Ml^ ^ {Ml + M|)2 - 4MIMI cos^ 2/3 
= Ml + M^ (5) 



Of course, to have a more accurate determination of all the Higgs boson masses and proper 
EWSB, it is important that the radiative corrections are included; see Refs. [2, 4] for reviews. 

In the non-universal Higgs mass (NUHM) model, the universality of the soft SUSY- 
breaking Higgs mass parameters at the GUT scale is relaxed, that is uihi 7^ '>tih2 7^ 'itT'O- 
This non-universality might, for instance, occur in GUT constructions where the fields Hi 
and H2 belong to different multiplets. This model is then equivalent to the "effective low 
energy" MSSM where in addition to the mSUGRA parameters, on has /j, and Ma as free 
input parameters. In this paper, we will parameterize this non-universality through two 
dimensionless parameters 5i and S2, which measure the relative deviation from the mSUGRA 
case at the GUT scale, as follows 

^^^(Mgut) = ml{MGVT){l + 5i) , m^,(MGUT) = ml{MGVT){l + ^2) (6) 
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from which one can clearly see that lower (larger) values of tuhj^ (^//a); corresponding to 
Si < {82 > 0), imply a lighter A boson. Thus, contrary to the mSUGRA case where the 
CP-odd A boson is very heavy for reasonable values of mo [except in the so-called focus 
point region where mo and tan/? are very large, sec Ref. [14] for instance], here, one can use 
the additional freedom and adjust the parameters 61^2 in such way that the obtained value 
of Ma_ is not too large. In particular, for high values of tan one can realize the scenario 
in which the Higgs sector of the model is in the intense coupling regime, i.e. Ma ~ 100-130 
GeV which leads to the relation Mjy fa Mh ~ Mh. 

Such a departure from universality leads to drastic consequences not only in the Higgs 
sector, but also in the chargino and neutralino sectors as the Higgs mass parameters m^j^ 
and enter in the determination of /i which enters as a basic input parameter in the mass 
matrices of these states. This would for instance alter the nature of the dark matter lightest 
neutralino with a possibly large impact on its cosmological relic density abundance and its 
detection rates, as will be summarized in the next subsections. 

2.2 Consequences on the LSP relic density 

In the MSSM there are four neutralinos, Xi with i — 1,2, 3, 4, among which the lightest one 
Xi is the hghtest SUSY particle (LSP) and the SUSY candidate for the cold dark matter 
in the universe .^J3eing a superposition of the bino, wino and higgsino fields, respectively 
denoted by B^, and H^, H2, the Xi neutralino texture can be written as 

xi = + Zi2W° + ZisH^ + ZuH^i . (7) 

where Zij are elements of the matrix which diagonalizes the 4x4 neutralino mass matrix. It 
is commonly defined that Xi is mostly gaugino-hke if P = + |^i2p > 0.9, higgsino-like 
if P < 0.1, and a mixed state otherwise. In the NUHM scenario, one can choose the values 
of niHi and mj^j at the GUT scale in such a way that the electrowcak symmetry breaking 
conditions lead to a low value for /i. As a consequence, the lightest neutralino is generally 
higgsino-like or a mixed bino-higgsino state. 

This has very important consequences on the relic density, which is inversely proportional 
to the LSP annihilation cross section. In Fig. 1, we display the Feynman diagrams of the 
main processes contributing to the annihilation cross section of the neutralinos to fermions 

(a), XiXi ~^ ffi s-nd to gauge or Higgs bosons (b), x?x? — >■ VV and with V = W, Z and 
$i = h,H,A and H^, together with the expressions of the relevant parts of the amplitudes 
[15]. From these diagrams and the LSP neutralino texture given by eq. (7), one can make 
the following remarks. 

• For LSP annihilation into light fermions [in general r-leptons and to a lesser extent 
b quarks], the diagram with t-channel sfermion exchange contributes only if Xi is 
bino-like [as higgsinos couple to fermions proportionally to their masses], while for a 
higgsino-like LSP, the coupling to the Z boson is large [being proportional to Z13 and 

Z14], enhancing thus the annihilation channel XiXi — ^ f f ■ 
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Figure 1: Feynman diagrams for LSP neutralino annihilation into a fermion pair (a) and 
into massive gauge bosons and Higgs bosons (b). The relevant parts of the amplitudes are 
shown explicitly. V and Z are the chargino and neutralino mixing matrices. 



The annihilation through Higgs boson exchange, XiXi 



A,h,H 



ff, contributes substan- 



tially only close to the pole of the Higgs boson exchanged in the s-channel and if the 
LSP is a mixed state in which case the coupling to the Higgs boson is large [both the 
Zii and Zi3{Zi4) matrix elements are large]; the dominant final state is bb for which 
the amplitude, oc nib tan (3, is enhanced for large values of tan /5. Note that the dom- 
inant component is due to the exchange of the pseudoscalar A boson as its s-wave 
contribution is not suppressed by the small velocity of the neutralinos, in contrast to 
the p-wave exchange of the CP-even h and H bosons. 

The annihilation into WW/ZZ bosons is efficient only for higgsino-like or mixed neu- 
tralinos when the elements Zy^iZx^ are large. There is also an annihilation diagram 
through ,s-channcl CP-cvcn /i, H bosons exchange but which gives only a small contri- 
bution as one is far from the Higgs boson poles and/or the Higgs couplings to neutralino 
and gauge bosons are suppressed; in addition, the p-wave contributions are suppressed 
by the small velocity of the LSPs. 

Since in the intense coupling regime all the Higgs particles are light, annihilation into 
a Higgs boson and a massive gauge boson, e.g. x'j'Xi ~^ or H^W^ , and two Higgs 
bosons, e.g. x^x? Ah or H^H~ , can occur if the neutralino LSP is heavy enough, 
m^o > 150 GeV. The dominant contribution is due to neutralino or chargino exchange; 
the Z{h, H) exchange diagrams give small contributions for Ah{AZ) annihilation as 
one is far from s-channel poles in this case. 

For higgsino-like LSPs, the co-annihilation processes with the lightest chargino and the 
next-to-lightest neutralino are very efficient because of the degeneracy m^± ~ m^o ~ 
m^o and lead to the right density only for very heavy LSPs [16]. There are additional 
diagrams which might contribute to the LSP cross section such as co-annihilation with 
the lightest sfermion [17] [in general fi] when it is almost degenerate with the LSP. 
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2.3 The impact on the detection of neutralinos 

Let us now discuss the expected signal of the intense coupling regime of the NUHM model 
in the two most promising search strategies for weakly interacting massive particle (WIMP) 
dark matter candidates: the search of the elastic scattering of ambient neutralinos off a 
nucleus in a laboratory detector through nuclear recoils, the "direct search" [18, 19], and the 
search for interesting decays products of WIMP annihilation into standard particles such as 
gamma rays, positrons, neutrinos etc., the so-called "indirect detection" [19, 20]. 

The strength of the direct detection signal is directly proportional to the neutralino- 
nucleon scattering cross section, a{xiN — >■ Xi^)- The matrix element for this scattering, 
mediated by squark and Z boson exchange as well as Higgs boson exchange diagrams [the 
crossed diagrams of those shown in Fig. la] receives both spin-dependent and spin-indepen- 
dent contributions. The former play a sub-dominant role in most direct search experiments, 
which employ fairly heavy nuclei. The spin-independent contribution in turn is usually 
dominated by Higgs boson exchange diagrams, where the Higgs bosons couple either directly 
to light {u, d, s) quarks in the nucleon, or couple to two gluons through a loop of heavy 
(c, b, t) quarks or squarks. Only scalar Higgs couplings to neutralinos contribute in the 
non-relativistic limit and therefore, in the absence of CP-violation in the Higgs sector of 
the NUHM model discussed here, one only needs to include the contributions of the two 
neutral CP-even h and H particles. The contribution of the CP-even Higgs boson which 
has enhanced couplings to down-type quarks for tan /3 1 is by far dominating. Thus, 
while in mSUGRA models only the contribution of the generally heavy H boson is relevant, 
in the intense coupling regime, both the h and H particles will play a role. 

Many direct detection experiments have been and arc presently carried out around the 
world. Recent collaborations such as CDMS [21] and EDELWEISS [22] have explored the 
regions of parameter space corresponding to a WIMP-nucleon cross section of a > 10^® pb. 
Future experiments will have improved sensitivity and, for example, GEDEON [23] will be 
able to explore a WIMP-nucleon cross section of cr > 3 x 10~^ pb. Similar cross sections will 
also be tested by the EDELWEISS II experiment while CDMS-Soudan [24] [an extension 
of the CDMS experiment in the Soudan mine], will be able to test a > 2 x 10^*^ pb [in our 
study, we will take the more conservative bound of 4 x 10~^ pb which corresponds to the 
present sensitivity of the experiment]. 

On the other hand, there are also promising methods for the indirect detection of WIMPs 
through the analysis of their annihilation into SM particles [with some of the contributing 
diagrams arc the same as those depicted in Fig. 1] and the search for their stable decay 
and fragmentation products, i.e. neutrinos, photons, protons, antiprotons, electrons and 
positrons [14]. While electrons and protons are undetectable in the sea of matter parti- 
cles in the universe, neutrinos, photons, positrons and anti-protons could be detected over 
the background due to ordinary particle interactions. An interesting possibility consists of 
detecting the gamma rays produced by these annihilations in the galactic halo. For this 
purpose, one uses atmospheric Cherenkov telescopes or space-based gamma-ray detectors. 
Planed experiments will reach significant sensitivity and, for instance, the GLAST telescope 
[25], which is scheduled for launch in 2007, will be able to detect a flux of gamma rays from 
dark matter particles of the order of $-y ~ 10~^° photons cm~^s~^. 
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The detection of the neutrahnos have been discussed at length in the mSUGRA scenario, 
see Ref. [8] for reviews. The rates can be increased in different ways when the structure of 
mSUGRA for the soft SUSY-breaking scalar terms is abandoned in the Higgs sector [18, 19, 
20]. In direct detection for instance, the neutralino-nuclcon clastic cross section, where the 
exchange of the CP-even Higgs bosons gives the dominant contribution as discussed above, 
can be considerably enhanced when decreasing the pseudoscalar Higgs mass and, thus, the 
heavy Higgs boson mass, through a particular choice of the non-universality parameters 
5i and 62- Concerning the prospects for the indirect detection of the neutralinos, we have 
seen that it is possible to enhance the annihilation channel involving the exchange of the 
CP-odd Higgs boson A by reducing its mass {Si < 0, ^2 > 0) and have a larger yield for 
the interesting decay products of the 6-quarks in the reaction XiXi ~* bb. In addition, by 
increasing the higgsino component of the lightest neutralino (^2 > 0), one would have a 
larger cross section for the annihilation processes XiXi ~^ WW/H^H~ for instance which 
allow to look for interesting decay products of the gauge and Higgs bosons such as photons. 

3. The analysis 

3.1 The determination of the spectrum and constraints 

In the present analysis, we use the Fortran code SuSpect [26] to solve the (two-loop) Renor- 
malization Group Equations for the gauge and Yukawa couplings and the soft SUSY-breaking 
parameters, to achieve proper electroweak symmetry breaking with the (two-loop) scalar 
potential, and to calculate the spectrum of the physical SUSY particles and Higgs bosons 
including radiative corrections. We follow the procedure outlined in Ref. [27], except that we 
allow for non-universality for the soft SUSY-breaking Higgs mass parameters as explained in 
the previous section. In addition to leading to a consistent electroweak symmetry breaking 
[i.e. obtaining a reasonable value for Ma and Mz as well as no charge or color breaking 
(CCB) minima] and to a phcnomcnologically viable spectrum [no tachyonic sparticlcs and 
Higgs bosons and and LSP that is not charged] , a given set of input parameters has to satisfy 
experimental constraints. The ones relevant for this study are as follows [28]. 

- The total cross section for the production of any pair of sparticles at the highest LEP 
energy (209 GeV) must be less than 20 fb. This leads for instance to an upper bound 
of 100 GeV on the chargino and slepton/ 3d generation squark masses. The masses 
of the gluino and first/second generation squarks should be larger than fa 200 GeV to 
cope with the Tevatron exclusion bounds. 

- Searches for neutral Higgs bosons at LEP impose a lower bound on Mh which, in 
the decoupling limit when ^ M^, is close to 114 GeV. Allowing for a theoretical 
uncertainty of p:^ 3 GeV, one requires > 111 GeV in this case. For small Ma values, 
the bound should be of the order of M^ ~ M^ > Mz- 

- Quantum corrections from superparticles to electroweak observables can be incorpo- 
rated through the p parameter which should obey the upper bound (^psusY < 2.2 • 10"^ 
at the 2fT level. However, it turns out that this constraint is always superseded by 
either the LEP Higgs search limit or by the CCB constraint. 



8 



- Recent measurements of the muon magnetic moment lead to a constraint on the SUSY 
contribution —5.7 • 10~^° < a^^^^ < 4.7 • 10~^ where, for the determination of the 
hadronic SM contributions, both data from e'^e" annihilation into hadronic final states 
and data from semileptonic r decays are used, leading to less than la deviation between 
the experimentally measured value and the theoretical prediction in the SM. 

- Allowing for experimental and theoretical errors, the branching ratio for radiative b 
decays should be 2.33 ■ 10^^ < B(6 — > 57) < 4.15 ■ 10~^. The branching ratio for the 
very rare decay into should be bounded by B{Bs — > n^n^) < 2.9 x 10~^ [29] 
[this observable does not yet constrain the parameter space of mSUGRA, but it it has 
been stressed recently [18] that it should be taken into account in the non-universal 
case]. Wc note however, that both the b ^ sj and Bg fi^fi" constraints have a 
different status from those discussed earlier, as a small amount of squark flavor mixing 
would make these observable compatible with the experimental values while having a 
negligible effect on the other (flavor conserving) observables, such as signals at coUiders. 

The theoretical constraints are implemented directly in the Fortran code SuSpect; the 
program implements also the experimental constraints mentioned above except for the rate 
B{Bs A*"*"/^") which we determine using the code microMEGAS [30]. 

In addition to these constraints from "collider experiments", we also require that the 
calculated Xi cosmological rehc density has to be in the the 99% confidence level WMAP 
narrow range, that is 

0.087 < l^^/i^ < 0.138 (8) 

where Q = p/pc with — 2 ■ 10~^^/i^g/cm^ is the "critical" mass density that yields a 
flat universe and the dimensionless parameter h the scaled Hubble constant describing the 
expansion of the universe. The result for the relic density has been obtained using the code 
microMEGAS [30]. 

We also take into account the most recent astrophysical bounds, in addition to the WMAP 
constraint. For the evaluation of the neutralino-nucleon cross section for direct detection and 
the gamma-ray fluxes relevant for the indirect detection, we use the latest released version 
of the program DarkSUSY [31]. We have included the relevant sensitivities and uncertainties 
in both detection methods. The measurement of the WlMP-nucleon cross section depends 
on the speed of the WIMP candidate in the neighborhood of the sun and it has been shown 
in Ref. [32] that the astrophysical uncertainties can affect the experimental sensitivities 
by a factor of two. Moreover, the modehzation of the 7r-nucleon can induce theoretical 
uncertainties in the calculation of the scattering WIMP-nucleon cross section, estimated 
to be a factor five at its maximum [33]; rather large QCD corrections and uncertainties 
are also present [34]. We will assume in this study an uncertainty factor of 2 coming from 
the astrophysical uncertainty on the WIMPs speed in our galaxy and a factor of 3 in the 
calculation of cr{xiN — > Xi^)- Concerning LSP indirect detection, most of the uncertainties 
are coming from the model used for the dark matter halo. For illustration, we will use in this 
paper a cuspy Navarro Frenk and White (NFW) profile [35], keeping in mind that a cuspier 
one [such as the one proposed by Moore in Ref. [36]] or a smoother (Isothermal) one, will 
give fluxes that are two orders of magnitude higher and lower, respectively. 
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3.2 The spectrum in the intense coupHng regime 

In the present study, we choose for iUustration two NUHM scenarios. In the first one, we 
fix tan/5 = 30, the common sfermion mass to mo = 2 TeV, the common gaugino mass to 
Ml/2 = 180 GeV and the triUnear coupUng in such a way that At = 3 TeV to maximize 
the mixing in the Higgs sector and thus the mass of the lightest h boson in order to pass 
easily the LEP2 constraint > 114 GeV in the decoupling regime when the theoretical 
uncertainty of ~ 3 GeV is included. We also fix the non-universality parameter to 5i = — 1, 
but vary the other parameter 62 in a narrow range to have a pseudoscalar A boson mass of 
the same order as the masses of the CP-even h and H bosons. In a second scenario, we 
choose a slightly lower value of tan (3, tan f3 = 20, and a larger value of the universal gaugino 
mass, Ml/2 = 580 GeV; all the other parameters are the same as in the previous case. 
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Figure 2: The neutral Higgs boson masses as a function of the parameter 62 in the two 
scenarios discussed in the text (top) and the cosmological relic density of the LSP neutralino 
(bottom). The pink vertical bands display the intense coupling regime while the brown 
horizontal lines show the WMAP range for the LSP relic density. 

In the upper part of Fig. 2, we display the masses of the three neutral Higgs bosons as a 
function of the parameter 62 in the two scenarios. The pink vertical bands correspond to the 
areas in which \MH,h — M^l < 10 GeV, that is, to the intense coupling regime which appears 
around 62 = 0.260 in the first case and 62 = 0.818 in the second case. These two values of 
62 lead to an almost mass degeneracy of the neutral Higgs particles Mh ~ Mh ~ Ma ~ 120 
GeV and define the points A and B corresponding to the initial values at the GUT scale : 
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niHi = and, respectively, mj^j = 1.26mo and mj^j = l-SlSmo- In the lower part of Fig. 2, 
we display the cosmological relic density of the LSP neutralino, again as a function of ^2, 
and compare it to the WMAP narrow band. For the values of 82 which correspond to the 
points A and B, we see that the WMAP constraint is obeyed. 

In the two scenarios, the output for the Higgs boson masses arc shown in Table 1 for the 
values of the parameters 82 which correspond to the selected points A and B. Wc also display 
the masses of the two lighter neutralinos m^o and m^o and the value of the parameter /x. As 
can be seen, in point A, the LSP x5 is rather light while the 11 parameter is very large: the 
LSP is thus almost bino-hke and the lighter chargino is almost mass degenerate with the 
wino-like while the heavier neutralinos and chargino are higgsino-like and have masses 
close to yU. In point B, the values of the gaugino and higgsino mass parameters are of the 
same order, M2 ~ yU, and all neutralinos and charginos are thus mixed states with masses 
that are quite close to each other. The LSP neutralino has a mass which is relatively large, 
exceeding the W/Z boson and all Higgs boson masses and even the top quark mass. 

In Table 1, we also display the result for the LSP relic density Q.^h?' and for the fractions of 
the various final states which contribute, B(x?Xi X) = a{xiXi X)/a{xiXi ~^ all)- As 
we have assumed heavy scalar fermions in both scenarios, mo = 2 TeV, the main annihilation 
channel for the LSP neutralinos, leading to a relic density compatible with the WMAP result, 
is the Higgs boson [mainly A and to a lesser extent h/H] exchange contributions in the case 
of point A. As can be seen in the table, in this case, the pseudoscalar A boson mass is close 
to twice the LSP neutralino mass and the Abb coupling is very strong for tan (3 — 30; this 
enhances the annihilation process XiXi ~* which reaches a branching fraction of 90 %. 
For point B with larger values of M1/2 and fi, higher values of the parameter 62 are needed to 
increase m^^j and keep a light pseudoscalar A boson, as can be seen from eq. (4). Since here, 
the A-pole is too far from the 2m^o threshold, the contribution of the annihilation process 
XiXi ~^ bb to the LSP rehc density is only at the level of 10%; the XiXi ~^ tt annihilation 
process, which has roughly the same probability as XiXi ~^ bb, occurs mainly through Z 
boson exchange. In fact, since in this scenario one has m^o ~ 250 GeV, all Higgs particles 
can occur as final states in the annihilation of the LSPs and, since the LSP is a bino-higgsino 
mixed state with strong couplings to the Higgs bosons, the fraction of the processes involving 
at least one final Higgs particle is larger than 70%. 

Finally, in Fig. 3, we show the result of a scan in the (M1/2, (^2) plane, starting with the 
basic inputs for tjiq = 2 TeV, tan/? = 30 (left) and 20 (right) and Si = —1; the position 
of the points A and B analyzed earlier are marked by a cross. In the white area, all the 
points fulfill the experimental and theoretical constraints discussed previously, except for the 
WMAP constraint which is obeyed only in the red narrow strips [between the strips, the relic 
density is either too high or too low, compared to the WMAP value] . The intense couphng 
regime, defined by the requirement iM^^h — Ma\ < 10 GeV, is represented by the thin black 
lines. In the (green) areas above these lines, the lightest h boson has a mass smaller than 
the LEP2 limit, Mh < 114 GeV. One notices that increasing the value of 62 for a given value 
of Ml/2, finally leads to a region in which M\ < : the mj^^ term is not sufficiently negative 
at the electroweak scale to compensate the mjj^ term in eq. (4) . 
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inputs 


M$ 






BR(x?X? ^ X) 


A 


tan/5 = 30 
62 = 0.260 
Ml/2 = 180 GeV 


Mh = 120.4 GeV 
Mh = 126.5 GeV 
Ma = 126.4 GeV 


m^o = 78 GeV 
m^o = 155 GeV 
/i = 846 GeV 


0.09 


bb : 90% 
rr : 10% 


B 


tan/3 = 20 
§2 = 0.818 
Ml/2 = 580 GeV 


Mh = 124 GeV 
Mh = 136 GeV 
Ma = 134 GeV 


m^o = 241 GeV 
m^o = 351 GeV 

A2 

/X = 362 GeV 


0.12 


bbitt) : 13%(12%) 
: 27%(1%) 
ZH{Zh) : 16%(5%) 
W^H^ : 22% 
PFiy : 2% 



Table 1: Sample spectra obtained in the non-universal scenarios A and B (with some of 
the parameters indicated and mo = 2 TeV), together with the value of the neutralino relic 
density and the fractions of the main subprocesses contributing to it. 




Figure 3: The allowed (M1/2, 62) parameter space in the NUHM model with tan/3 = 30 
(left) and tan/3 = 20 (right), /i > 0, A = 3 TeV, mo = 2 TeV and 61 = -1. We indicate 
in the plot the points A and B defined in the text, the WMAP range, the area in which 
Mh < 114 GeV, and the intense coupling regime (black solid lines). 



3.3 Prospects for dark matter detection 

Let us now turn to the prospects for detecting the LSP dark matter candidate in these 
two scenarios. In direct detection, as a consequence of the hght Higgs bosons with strongly 
enhanced couplings to down-type quarks that are present, the scalar or spin-independent 
LSP-nucleon scattering cross section can be extremely large as the process is mainly me- 
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diated by the exchange of the two CP-even Higgs bosons h and H. As can be seen in 
the upper part of Fig. 4, where the expected Xi-^ ~^ Xi^ scattering cross section and the 
CDMS-Soudan sensitivity in which the various uncertainties discussed in section 2.3 are in- 
corporated, the intense-couphng regime is within the reach of the CDMS-Soudan detector. 
For instance, in point A with 62 = 0.260, one has a{xiN — > Xi^) = 1-2 x 10"''' pb which is 
only shghtly lower than the present CDMS sensitivity of 4 x 10~^ pb and much larger than 
the expected sensitivity of 2 x 10~^ pb. In the second scenario, one obtains an even larger 
cross section despite of the larger value of the LSP mass and the smaller value of tan jS. For 
instance, in point B with 62 = 0.818, one obtains a{xiN) = 10^^ pb. This is due to the fact 
that here, the LSP is a mixed bino-higgsino state with larger couplings to the Higgs bosons 
[in the first scenario, the larger value of tan/3 and the smaller LSP mass were partly com- 
pensated by the smaller Higgs couplings to the almost bino-like neutralino]. Note that such 
high values for the LSP-nucleon scattering cross section have also been obtained in Ref. [18] 
in the case of light pseudoscalar Higgs bosons not pertaining to the intense-coupling regime. 




-10 




0.825 



Figure 4: Upper curves: the rate for the LSP-nucleon cross section including the various er- 
rors (hatched bands) compared to the expected sensitivity of the CDMS-Soudan experiment 
(horizontal blue bands) as a function of the non-universality parameter 62- Lower curves: 
the obtained gamma ray flux in photons cm^^s^^ using the NFW profile as a function of 
S2 (dashed curves) and the expected sensitivity of the GLAST experiment (dotted curves). 
The left (right) handed plots are for the scenarios discussed previously with M1/2 = 180 GeV 
(580 GeV) and the other parameters as described in the text. 
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Concerning indirect LSP detection, we have calculated the rate in the framework of the 
NUHM model, assuming the NFW profile [35] within a sohd angle of AQ = 10~^ sr and 
found in the two scenarios fluxes around or larger than ~ lO^^^cm^^s"^, which is within the 
reach of a detector like GLAST. This is shown in the lower parts of Fig. 4, where wc display 
the obtained gamma ray flux as a function of 82 and compare it to the expected sensitivity 
of the experiment. To be more precise, we find a flux of 5.6 x 10~^°cm~^s~^ for 

point A and _ g 4 ^ 10~^^cm~^s~^, for point B. In the first case, the obtained fiux is 

much larger than in the second case because of the stronger enhancement of the annihilation 
channel through 74-boson exchange and the smaller value of the LSP neutralino mass. 

Thus, one can conclude that the two scenarios that we have chosen to illustrate the Higgs 
intense coupling regime in the NUHM model are in a sense orthogonal, as point A will be 
more easily probed in indirect detection experiments whereas point B can be best probed in 
direct LSP detection. 

4 The mSUGRA case 

For completeness, let us briefiy discuss the intense coupling regime in the mSUGRA frame- 
work with universal Higgs mass parameters. As mentioned earlier, the pseudoscalar Higgs 
boson mass and the // parameter tend to be quite large in this case, except in the focus point 
region when both the common scalar mass tjiq and tan/3 are large. This, however, occurs 
only in a tiny part of the parameter space near the region where electroweak symmetry 
breaking is not achieved [in general because either M\ < or /i,^ < or both]. The Higgs 
intense coupling regime occurs thus in an even smaller region of the parameter space and 
one has to tune severely some of the input parameters to realize it. 

For some of the input parameters already used in the NUHM model, mo = 2 TeV and 
At = "i TeV, and assuming fi > and M1/2 = 350 GeV for instance, one obtains using 
the program SuSpect, the intense coupling regime only for values of tan/? very close to 
tan/5 = 53.87 as is exemplified in Fig. 5, where the neutral Higgs boson masses are shown 
in this mSUGRA scenario as a function of tan /3. For the particular value tan /3 = 53.87, 
one has a pseudoscalar Higgs boson mass of Ma — 122 GeV leading to Mh ~ 117 GeV 
and Mh — 124 GeV for the masses of the CP-even Higgs particles. For the chosen inputs, 
one also obtains a relatively large value for the higgsino mass parameter, fi ~ 440 GeV, 
with a lightest neutralino mass of m^o ~ 150 GeV. The LSP is thus almost bino-like with 
only a small higgsino component. However, this component is large enough to lead to a 
substantial annihilation rate of the LSP neutralino into Higgs bosons. Indeed, one obtains 
for instance B{XiXi ~^ hA + HA) and B{xiXi ~^ HZ + hZ) of, respectively, 12% and 
8% the dominant channel being annihilation into hb and t^t~ final states which has an 
overall fraction close to 70%. The obtained value for the LSP cosmological relic density in 
this scenario is Vt^h? ~ 0.1 which is within the WMAP range. Because of the very high 
value of tan/3, the couphngs of the rather fight CP-even neutral Higgs bosons h and H 
to hb states are extremely strong, leading to a large LSP-nucleon scattering cross section, 
(^{Xi^ ~^ Xi^) ~ 2 X 10~^ pb, which can be easily probed by the CDMS experiment. The 
indirect detection rate of the lightest neutralino is also significant as one obtains a gamma- 
ray flux of the order of ^J^fw ^ 2 x lO^^^cm^^s"^ which can be observed by GLAST. 
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Figure 5: The same as in Figs. 2 and 4 but for the mSUGRA model with universal Higgs 
mass parameters and as a function of tan /3 with basic inputs mo = 2 TeV, At = 3 TeV, 
Ml/2 = 350 GeV and > 0. 

However, there is a potential problem with this particular scenario. For such a large value 
of tan/? and a low pseudoscalar Higgs boson mass Ma, the decay rate B{Bs — > oc 
tan® is huge [as already noticed in Ref. [37] for instance] and one obtains for the point 

introduced above B{Bs ~ 4 x 10~^ to be compared with the Tevatron experimental 

bound of B{Bs —>■ fj^^fi'^) ~ 2.9 x 10^^ [29]. As discussed in section 3.1, the two-orders of 
magnitude discrepancy between these two values will need a rather strong flavor mixing to 
be resolved. [Note that we have a similar problem for the point A of the NUHM model, but 
the difference between the obtained Bg — > /x'''/^" rate and the experimental bound is much 
smaller and can be attributed to theoretical uncertainties.] 

Nevertheless, if this problem can be circumvented, one would be very close to definitely 
test this scenario. Indeed, the cross section for the associated production of the neutral 
Higgs bosons in association with bb pairs, pp —>■ gg/qq —>■ bb + A,h,H, will be extremely 
large at the Tevatron for the considered values of tan/3 and M4. With the luminosity that 
is presently collected at the Tevatron, values tan/3 ~ 50 are already ruled out for Ma ~ 100 
GeV [38] and the 2007 projection for the increase of the Tevatron luminosity would allow to 
observe or rule out a Higgs boson signal in this process for pseudoscalar Higgs bosons masses 
of Ma ~ 120 GeV [39]. 
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5. Conclusions 



In this paper, we have investigated the possibihty of reahzing the intense couphng regime, 
in which all the Higgs bosons of the MSSM have comparable masses and the parameter 
tan/3 is large, in the framework of the constrained MSSM. We have pointed out that this 
interesting regime is possible within the mSUGRA model, but when all the scalar fermion 
and Higgs soft SUSY-breaking mass terms are unified at the GUT scale, this happens in a 
very limited range of the parameter space: large values of tan f3, > 50, are needed which 
might lead to a severe conflict with the measured value of the Bg fi^fi'^ branching rate and 
potentially, with the search of the MSSM Higgs bosons at the Tevatron. In contrast, relaxing 
the universality of the SUSY-breaking scalar masses of the two MSSM Higgs fields at the 
high scale, niHi 7^ ^j?2 ''^Oi leads to more freedom which allows to reach this regime in a 
much larger area of the mSUGRA parameter space. We have provided examples of scenarios 
in which the intense coupling regime is realized in the NUHM model and where, not only 
all known theoretical requirements and collider experiment constraints are fulfilled, but also 
where the generated rehc density of the lightest neutrahno of the MSSM is compatible with 
the WMAP measurement of the amount of cold dark matter in the universe. 

These scenarios have a very interesting phenomenology. They lead by definition to light 
MSSM Higgs bosons which can be observed at the next generation of colliders such as 
the LHC and the ILC [and also at the Tevatron if the value of tan/3 is very large]. In 
addition, because of the large tan (5 and low Ma values which are implied by these models, 
the rates for direct and indirect detection of the LSP cold dark matter neutrahno are large 
enough for this particle to be observed in near future experiments such as CDMS-Soudan 
and GLAST. Furthermore, since at least the spectrum of spin-1/2 superparticles needed to 
generate the right cosmological density for the LSP is relatively light, the heavier charginos 
and neutralinos [and possibly the gluinos] could also be produced at the next generation of 
collider experiments. This would allow to precisely measure the couplings of these particles, 
in particular the couplings to the Higgs bosons in possible cascades [40], and determine the 
parts of the MSSM Lagrangian which enter in the expressions of the relic density and the 
detection rates. Thus, using the future collider data, one would be able to check that the 
LSP neutrahno makes indeed the entire dark matter of the universe and that the detection 
rates observed in astrophysical experiments are indeed those predicted by the model. In this 
respect, these scenarios, if realized in Nature, would highlight the complementarity between 
collider physics searches and precision measurements and astroparticle physics searches. 
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